r] —>■ 7r°77 decay within a chiral unitary approach revisited 



OO 
O 
O 



Oh' 
I 

Oh' 

CD 



> 

cn 
cn 

O 

00 

o 



X 



E. Oset\ J. R. Pelaez^ and L. Roca^ 
^ Departamento de Fisica Tedrica and IFIC, Centra Mixta Universidad de Valencia-CSIC, 
Institutas de Investigacion de Paterna, Aptdo. 22085, 46071 Valencia, Spain 
^ Departamento de Ftsica Teorica 11, Universidad Complutense. 28040 Madrid, Spain. 
^ Departamento de Fisica, Universidad de Murcia, E-30071, Murcia, Spain 
(Dated: February 2, 2008) 

In view of the recent experimental developments on the experimental side in the r; 7r°77 decay, 
and the fact that the Particle Data Group in the on line edition of 2007 reports sizable changes of 
the radiative decay widths of vector mesons used as input in the theoretical calculations of [l|], a 
reevaluation of the decay width of the 77 in this channel has been done, reducing its uncertainty by 
almost a factor of two. The new input of the PDG is used and invariant mass distributions and 
total widths are compared with the most recent results from AGS, MAMI and preliminary ones of 
KLOE. The agreement of the theory with the AGS and MAMI data is very good, both for the total 
rates as well as for the invariant mass distributions of the two photons. 

PACS numbers: 13.40Hq, 12.39Fe 



I. INTRODUCTION 

The 77 — > 7r°77 reaction has been quite controversial 
given the large discrepancies between different theoret- 
ical approaches trying to match the scarce experimen- 
tal data. For a long time the standard experimental re- 
sults have been those of early experiments giving 
r = 0.84 ± O.lSeV. More recent experiments with the 
Crystal Ball detector at AGS Q reduced this value to 
r = 0.45 ± 0.12 eV. A new reanalysis of AGS data gives 
r = 0.285 ± 0.031 ± 0.049 eV |5i] and a more recent anal- 
ysis with the Crystal Ball at MAMI provides the rate 
r = 0.290 ± 0.059 ± 0.022 eV d. At the same time the 
last two experiments have provided the much awaited 
invariant mass distribution for the two photons, which 
was thought to provide valuable information concerning 
the theoretical interpretation. Some preliminary results 
from KLOE at Frascati @ are also available with values 
around T = 0.109 ± 0.035 ± 0.018 eV. 

The theoretical models show also a similar dispersion of 
the results, from large values obtained using models with 
quark box diagrams to much smaller ones, obtained 
mostly using ideas of chiral perturbation theory (ChPT) , 
which are quoted in jlj. 

The r] —> ir'^jj reaction has been traditionally consid- 
ered to be a border line problem to test chiral pertur- 
bation theory (ChPT). The reason is that the tree level 
amplitudes, both at 0{p^) and 0(j5^), vanish. The first 
non- vanishing contribution comes at O(p^), either from 
loops involving kaons, largely suppressed due to the kaon 
masses, or from pion loops, again suppressed since they 
violate G parity and are thus proportional to m,j — 
p. The first sizable contribution comes at O(p^) but the 
coefficients involved are not precisely determined and one 
must recur to models. In this sense, either Vector Meson 
Dominance (VMD) [l^, [ll[ , the Nambu- Jona-Lasinio 
model (NJL) 1131 . or the extended Nambu- Jona-Lasinio 
model (ENJL) [ll,[l3, have been used to determine these 
coefficients. However, the use of tree level VMD to ob- 



tain the 0{p^) chiral coefficients by expanding the vector 
meson propagators, leads [1] to results about a factor of 
two smaller than the "all order" VMD term when one 
keeps the full vector meson propagator. The lesson ob- 
tained from these studies is that ChPT can be used as a 
guiding principle but the strict chiral counting has to be 
abandoned since the 0{p^) and higher orders involved in 
the full ( "all order" ) VMD results are larger than those of 
0{p*). Also these calculations had several sources of un- 
certainty, one of the most important was the contribution 
of the ao(980) resonance, for which not even the sign was 
known. Thus, one is lead to rely directly on mechanisms 
for the reaction, leaving apart the strict chiral counting. 

The theoretical situation improved significantly with 
the thorough revision of the problem in [1| , where the dif- 
ferent sources of uncertainty were studied and the ao (980) 
contribution was reliably include d by using the unitary 
extensions of ChPT [lE HI, |l3, Within this chi- 

ral unitary approach for the interaction of pseudoscalar 
mesons the ao(980), as well as the /o(980) or the cr(600) 
resonances, are dynamically generated by using as input 
the lowest order chiral Lagrangians (Toj and resuming the 
multiple scattering series by means of the Bethe Salpeter 
equation. Another source of corrections in 1] was the use 
of the newest data for radiative decay of vector mesons 
of the PDG 2002 Q. It was noted in [1] that the rates 
had significantly changed from previous editions of the 
PDG, to the point that the rj 77^77 widths calculated 
in [l3| would have changed by about a factor of two 
should one have used the new data for radiative decay 
of vector mesons of the PDG 2002 instead of the former 
ones. Another improvement in [l| was the unitarization 
of the pair of mesons of the VMD terms beyond the tree 
level. Furthermore, to have a better control on the reac- 
tion, the consistency of the model with the related reac- 
tion 77 7r°77 was established. Finally, in [ij a thorough 
analysis of the theoretical errors was done by considering 
all sources of uncertainty and making a MonteCarlo sam- 
ple of results obtained with random values of the input 
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within the uncertainties. 

The final resuh obtained in [l[ was 

r = 0.42 ±0.14 eV, 



(1) 



which is stiU in agreement with the present experimen- 
tal results within uncertainties. Nevertheless, five years 
after the publication of these results some novelties have 
appeared that call for a revision of the problem. In- 
deed, once again the data for the radiative decay of vector 
mesons of the "on line" PDG 2007 20] have significantly 
changed with respect to the data of the PDG 2002 used 
in ref. [l| . The correction due to these changes is impor- 
tant and it produces about a 25% decrease in the central 
value of the result of Eq. ([1]). At the same time, the 
theoretical uncertainty is reduced by almost a factor of 
two. On the other hand, the new experimental results 
regarding the two photon invariant mass distribution 
provide an extra challenge for the theoretical models. 

In view of this, it has become necessary to update the 
work of [l| to account for the newest experimental results 
of the PDG 2007 [l^] and to compare with the most 
recent experimental data of the 77 — *■ 77*^77 decay. The 
model used here is, hence, the same as the one of [1] 
and the only changes are the use as input of the new 
vector mesons radiative widths. Thus, we refrain from 
providing detailed explanations on the model and in this 
brief report we just concentrate on the changes. 



II. VMD CONTRIBUTION 

Following we consider the VMD mechanism of Fig. [T] 
which can be easily derived from the VMD Lagrangians 




for the different radiative decays, together with the the- 
oretical results (using Gy = 69MeV and / = 93MeV) 
and experimental [3, [131 branching ratios. In Table IJ 
we quote the results of the PDG version of 2002, which 
were used as input in the evaluation of the results in [l| , 
together with the new results of the PDG 2007 on-line 
edition j2^ which are used in the present paper. 

The agreement of the theoretical results with the data 
is fair but they can be improved by incorporating SU (3) 
breaking mechanisms ^2^. For that purpose, we normal- 
ize here the Ci couplings so that the branching ratios in 
Table U agree with experiment. 

Once the VP^ couplings have been fixed in this way, 
we can use them in the VMD amplitude corresponding 
to the diagram of Fig. [U for what we follow the details 
of [1]. Next we briefiy describe the other mechanisms 
considered in 



III. OTHER MECHANISMS 

In other mechanisms were considered which are not 
affected by the modifications of the previous section. We 
refresh them graphically and forward the reader to [1] for 
details. 

In Fig. [5] we show the diagrams that go through kaon 
loops. These diagrams, with the unitarization of the 
meson-meson interaction depicted in Fig. [3l were shown 
in [23] to be mostly responsible for the strength of the 
77 7r°77 reaction in the region of the ao(980) reso- 
nance. It was also shown in [l[ that the consideration of 
the mechanisms of Fig. [T] and Fig. d] improved the agree- 
ment with the data at low 77*^77 invariant masses. 



FIG. 1: Diagrams for the VMD mechanism, 
involving VVP and couplings [2l[ 

Cvvp = ^e^''""'{d^V,dc.VpP), Cv-, = -^fegA^{QV^' 

(2) 

where and P are standard SU{2>) matrices for the 
vector mesons and pseudoscalar mesons respectively [l|. 

In Eq. © G = J^, 5 = m and / = 93MeV, 

with Gy the coupling of p to tttt in the normalization of 
[2^ . From Eq. ^ one can obtain the radiative decay 
widths for V P7, which are given by 



FIG. 2: Diagrams for the chiral loop contribution 
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(3) 



where k is the photon momentum for the vector meson at 
rest and Ci are SU (3) coefficients that we give in Table U 



FIG. 3: Resummation for 77 ix^r). 

The vector meson exchange diagrams of Fig. [T] were 
unitarized in [ij by including the resummation of dia- 
grams of Fig. [21 producing the diagram depicted in Fig. 31 
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TABLE I: SU(3) d coefficients together with theoretical and experimental branching ratios for different vector meson decay 
processes. 
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where the thick dot represents the fuU meson-meson uni- 
tarized amphtude. Note that these mechanisms are also 
affected by the renormahzation of the VVP vertices dis- 
cussed in the previous section. 
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FIG. 4: Loop diagrams for VMD terms. The diagrams with 
the two crossed photons are not depicted but are also included 
in the calculations. 

Finally, a small term related to the three meson axial 
anomaly, and shown diagrammatically in FigJSl was also 
taken in the calculation since, as noted in [9|, although 
small by itself gives a non-negligible contribution upon 
interference with the other terms. 




FIG. 5: Diagrams with two anomalous 7 3M vertices. 



IV. RESULTS 

By considering all the modifications discussed in sec- 
tion ini the integrated width that we obtain is 



r = 0.33 ±0.08 eV 



(4) 



which should be compared to the result of 1] of F = 
0.42 ± 0.14 eV. The new result compares favorably with 
the most recent results of Cristal Ball at AGS F = 0.285± 



0.031±0.049 eV and MAMI F = 0.290±0.059±0.022 eV 
Q. However, all these decay widths are much larger than 
the preliminary results of KLOE at Frascati F = 0.109 ± 
0.035 ±0.018 eV @. 

The mass distribution of the two photons provides ex- 
tra information which was claimed to be relevant to fur- 
ther test theoretical models. In [l| the differential cross 
section dV / dM^^ was given. We present here the updated 
results in Fig. [6l where the contribution of the different 
mechanisms is shown. The new experiments reported in 
provide measurements of dT/dM^^ which can be con- 
trasted with theoretical predictions. 




FIG. 6: Contributions to the two photon invariant mass distri- 
bution. From bottom to top, short dashed line: chiral loops; 
long dashed line: only tree level VMD; dashed-dotted line: co- 
herent sum of the previous mechanisms; double dashed-dotted 
line: idem but adding the resummed VMD loops; continuous 
line: idem but adding the anomalous terms of Fig. [5] which is 
the full model presented in this work (we are also showing as a 
dotted line the full model but substituting the full tx+if- ,r]TvO 
amplitude by its lowest order). 

Note that in the experiments of Q the magnitude 
dr/dM^^ is given, while in |l[ and in Fig. MdT / dMj^ is 
evaluated. Although these distributions are equivalent, 
in practice the first one is more useful to study the spec- 
trum at low invariant masses since it provides extra in- 
formation not given by the second one. Indeed, dV /dM^^ 
is zero at the threshold of the 77 phase space. However, 
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dV/dM"^^ contains an extra l/2M-y^ factor and leads to 
a finite value at zero 77 invariant mass. This finite value 
and the shape of the distribution close to threshold offer 
an extra test to the theory that would be missed had we 
simply taken dV/dM^ry for comparison. This of course 
implies that the measurements can be done with good 
precision at the threshold. On the other hand, for the 
high mass region of the spectrum the dT/dMjj distri- 
bution is more suited to reveal the effects of different 
theoretical mechanisms, as we have shown in Fig. [S) 

In Figs. 7(a) and 7(b) we compare the theoretical re- 
sults that we obtain with the distributions obtained at 
MAMI and AGS. The agreement is good, in shape and 
size, and the theory provides indeed a finite value at 
threshold compatible with experiment, which has nev- 
ertheless large errors. It is interesting to see that the 
AGS data show clearly an increase of the distribution at 
low invariant masses which is a feature of the theoretical 
results. The data of MAMI, however, have too large er- 
rors at threshold and does not allow one to see this trend 
of the results. At large values of the invariant mass the 
agreement of the theory with the MAMI data is better 
than with the AGS data. 

In order to offer a different perspective of the com- 
parison of the results at the higher mass region of the 
distribution, we show in Figs. 7(c) and [7(d)] our final re- 
sults for dT/dM-y-y compared to the data of ^] properly 
transformed to these variables. 



V. SUMMARY 

In summary, we have witnessed an important exper- 
imental advance in the recent years on the rj 77^77 



decay. The parallel advances in theory reflected by the 
work of [2] have allowed a detailed comparison of results 
which has given a good agreement both for the total rate 
as well as for the invariant mass distributions with the 
most recent finished results. The discrepancy with the 
preliminary data of Frascati is worrisome, but we should 
wait till these data are firm before elaborating further on 
the discrepancies. 
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